Far-UV FUSE spectroscopy of the O VI resonance doublet in 

Sand 2 (WO) 

Paul A. Crowther 1 , A. W. Fullerton 2,3 , D. J. Hillier 4 , K. Brownsberger 5 , L. Dessart 1,6 
A. J. Willis 1 , 0. De Marco 1 , M. J. Barlow 1 , J. B. Hutchings 7 , D. L. Massa 8 , 

D. C. Morton 7 , G. Sonneborn 9 

Received 17 March 2000; accepted XX May 2000 
Submitted to Ap.J. Letters 



department of Physics & Astronomy, University College London, Gower Street, London 
WC1E 6BT, UK 

2 Dept. of Physics & Astronomy, University of Victoria, P.O. Box 3055, Victoria, BC, 
V8W 3P6, Canada. 

3 Department of Physics & Astronomy, Johns Hopkins University, 3400 North Charles St., 
Baltimore, MD 21218 

4 Department of Physics & Astronomy, University of Pittsburgh, 3941 O'Hara Street, PA 
15260 

5 Center for Astrophysics & Space Astronomy, Campus Box 389, Boulder, CO 80309 
6 Present Address: Departement de Physique, Universite Laval and Observatoire du Mont 

Megantic, Quebec, QC, Canada, G1K 7P4 

7 D.A.O., Herzberg Institute of Astrophysics, National Research Council of Canada, 5071 

W. Saanich Road, Victoria BC, Canada, V8X 4M6 

8 Raytheon ITSS, NASA's Goddard Space Flight Center, Code 631, Greenbelt, MD 20771 
9 Laboratory for Astronomy and Solar Physics, NASA/GSFC, Code 681, Greenbelt, MD 

20771 



- 2 - 



ABSTRACT 

We present Far-Ultraviolet Spectroscopic Explorer (FUSE) spectroscopy of 
Sand 2, a LMC WO-type Wolf-Rayet star, revealing the VI resonance P Cygni 
doublet at 1032-38A. These data are combined with HST/FOS ultraviolet and Mt 
Stromlo 2.3m optical spectroscopy, and analysed using a spherical, non-LTE, line- 
blanketed code. Our study reveals exceptional stellar parameters: T* ~150,000K, 
^00=4100 kms" 1 , log(L/L )=5.3, and M = 1 x lO^Moyr" 1 , if we adopt a 
volume filling factor of 10%. Elemental abundances of C/He~ 0.7 ± 0.2 and 
0/He~ 0.151^5 by number qualitatively support previous recombination line 
studies. We confirm that Sand 2 is more chemically enriched in carbon than 
LMC WC stars, and is expected to undergo a supernova explosion within the 
next 5x 10 4 yr. 

Subject headings: stars: Wolf-Rayet - stars: evolution - stars: individual 
(Sand 2) 

1. Introduction 

Wolf-Rayet (WR) stars provide keys to our understanding of massive stellar evolution, 
nucleosynthesis processes and chemical enrichment of the ISM. Of these, the oxygen- 
sequence (WO) introduced by Barlow & Hummer (1982) is by far the rarest. Their 
high excitation oxygen emission lines are widely interpreted as revealing the late core 
helium-burning or possibly carbon-burning stage (Smith & Maeder 1991), of importance for 
constraining the controversial 12 C(a,7) 16 reaction rate. 

In spite of searches in all Local Group galaxies, only six massive WO stars are known 
to date, namely Sand 1 (Skl88) in the SMC, Sand 2 (BAT99-123) in the LMC, Sand 4 
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(WR102), Sand 5 (WR142) and MS4 (WR30a) in our Galaxy, and DR1 in IC1613. Since 
Ovi 38 11-34 A is a primary WO classification diagnostic, with an equivalent width of up to 
1700A (Kingsburgh, Barlow & Storey (1995, hereafter KBS), observations of O VI 1032-38A 
are keenly sought. However, its location in the far-UV has ruled out such observations 
to date. This situation has changed following the successful launch of the Far- Ultraviolet 
Spectroscopic Explorer (FUSE, Moos et al. 2000), which permits routine high dispersion 
far-UV spectroscopy of massive stars in the Magellanic Clouds. 

In this Letter, we analyse FUSE spectroscopy of Sand 2 (Sanduleak 1971), alias Sk— 68° 
145 = Brey 93 = BAT99-123 (Breysacher, Azzopardi & Testor 1999), together with Hubble 
Space Telescope (HST) and ground-based datasets. 

2. Observations 

Previously unpublished far-UV, UV and optical/near-IR spectroscopy of Sand 2 have 
been obtained with FUSE, HST and the Mt Stromlo and Siding Spring Observatory 
(MSSSO) 2.3m, respectively. 

2.1. Far-UV spectroscopy 

Sand 2 was observed by FUSE as part of the Early Release Observation programme 
X018 on 1999 Oct 31. A 8134 sec exposure of Sand 2 with the 30" x 30" (LWRS) aperture 
provided data at R ~12,000 with the two Lithium Fluoride (LiF) channels, covering 
AA979-1187A, obtained in time-tag (TTAG) mode. At this epoch, the two Silicon Carbide 
(SiC) channels, covering AA905-1104A, were badly aligned so that SiC data were of poor 
quality. 
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Sand 2 data were processed through the pipeline data reduction, CALFUSE (version 
1.6.8), and are shown in Fig. 1. The pipeline extracted the ID spectrum, removed the 
background, and corrected for grating wobble and detector drifts. No corrections for 
astigmatism or flat fielding have been applied. From Fig. 1, the far-UV continuum flux 
of Sand 2 is low 5 x 1CT 14 ergs -1 cm -2 A -1 ), with two principal stellar features, Ovi 
A1032-38 and Civ A1168. The latter is blended with Cm A1175, while Civ A1108 and 
O VI A1125 are present, though weak. The FUSE spectrum is affected by a multitude of 
interstellar absorption features, principally Hi and H 2 , plus airglow emission features due 
to Ly(3, [Ol] and [Nl]. 

2.2. Near-UV spectroscopy 

Sand 2 was observed with the HST Faint Object Spectrograph (FOS) instrument 
during 1996 March (PI: D.J. Hillier, Program ID 5460). Exposures totalling 5980, 3780 
and 1090 sec were obtained with the G130H, G190H and G270H gratings, respectively. 
Many important spectral features are identified in the HST/FOS dataset (see also KBS), 
principally Civ A1550, Cm A2297, Hen A1640, Oiv A1340, Ov A1371, A2781-87. 
Following Prinja, Barlow & Howarth (1990) we find foo=4100 kms -1 from Civ A1548-51 
(KBS derived 4500 kms" 1 from IUE spectroscopy). 

2.3. Optical spectroscopy 

We have used the Double Beam Spectrograph (DBS) at the 2.3m MSSSO telescope to 
observe Sand 2 on 1997 Dec 24-27. Use of a dichroic and the 300B, 600B and 316R gratings 
permitted simultaneous spectroscopy covering 3620-6085A and 6410-8770A for 1200 sec, 
plus 3240-4480A and 8640-11010A for 2500 sec. A 2" slit and 1752x532 pixel SiTE CCD's 



- 5 - 



provided a 2 pixel spectral resolution of ~5A. A standard data reduction was carried out, 
including absolute flux calibration, using wide slit spectrophotometry of HD 60753 (B3 IV) 
and \x Col (09 V), plus atmospheric correction using HR2221 (B8V). Convolving our 
dataset with Johnson broad-band filter profiles reveals V=15.1 and B-V=+0.54 mag. Since 
these are contaminated by emission lines, we have convolved our data with Smith ubvr 
narrow-band filters to reveal v—16.1, u — b = —0.19, b — v = —0.07 and v — r — —0.14 mag. 

Our optical dataset confirms the spectral morphology previously presented and 
discussed by KBS, with exceptionally strong Oiv A3400, Oiv A3818, Civ+Hell A4660 and 
Civ A5801. Depending on the WO selection criteria, its spectral type is W03 (Crowther, 
De Marco & Barlow 1998) or W04 (KBS). 

3. Quantitative analysis of Sand 2 

Wolf-Rayet winds are so dense that non-LTE effects, spherical geometry and an 
expanding atmosphere are minimum assumptions. 

3.1. Analysis technique 

We employ the code of Hillier & Miller (1998), cmfgen, which iteratively solves the 
transfer equation in the co-moving frame subject to statistical and radiative equilibrium in 
an expanding, spherically symmetric and steady-state atmosphere. These models account 
for clumping, via a volume filling factor, /, and line blanketing, both of which have a 
significant effect on the physical properties of WC stars (e.g. Hillier & Miller 1999). 
Through the use of 'super-levels', extremely complex atoms can be included. For the 
present application, a total of 3,552 levels (combined into 796 super- levels) , 60 depth points 
and 63,622 spectral lines of Hel-ll, Cll-iv, O il-vi, Nell-lV, Siiv, SlV-vi, Arlll-v, and 



- 6- 



Fe IV- vill are considered simultaneously (see Dessart et al. 2000 for the source of atomic 
data used). 

We adopt a form for the velocity law (Eqn 8 from Hillier & Miller 1999) such that 
two exponents are considered (/3i=l, (32=50, t> cx t=2900 kms -1 , f tU rb=100 kms" 1 ), with the 
result that acceleration is modest at small radii, but continues to large distance (0.9-Uoo is 
reached at 100i?*). 

Our spectroscopic analysis derives M/y 7 /, rather than M and / individually since line 
blending is severe in Sand 2. A series of models were calculated in which stellar parameters, 
T*, log(L/L ), M/v7) C/He and O/He, were adjusted until the observed line strengths 
and spectroscopic fluxes were reproduced. We adopt 0AZ & abundances for Ne, Si, S, Ar 
and Fe. The distance to the LMC was assumed to be 51.2 kpc (Panagia et al. 1991). 

The wind ionization balance is ideally deduced using isolated spectral lines from 
adjacent ionization stages of carbon (Cm A2297, Civ A5801-12) or oxygen (Oiv A3404-14, 
O V A3144, O VI A5290). In practice this was difficult to achieve for Sand 2 because of the 
severe line blending. 

3.2. Stellar parameters 

Our initial parameter study of Sand 2 revealed a fairly similar spectral appearance 
spanning 120kK< T, <170kK, with log (L/L )=5.28 and log (M/A^yr" 1 ) = -4.94 fixed. 
The principal differences between these models are that (i) the observed strength of the 
Ovi 3811-34A doublet, plus lines in the red such as Civ 7700A favour a high T*; (ii) the 
weakness of the Ovi 1032-38A doublet, as revealed by FUSE, favours a lower T*. Since 
other parameters, in particular abundances, are largely unaffected by these discrepancies, 
we shall adopt T*=150kK (i.e. i?*=O.65-R ). Note that higher luminosity models do 
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produce significant effects, such as a dramatic weakening of Civ 580 1-12 A emission. 

Fig. 2 compares our synthetic spectrum with the observed far-UV, UV and optical 
spectroscopy of Sand 2. Our model is reddened by E(B — U)=0.08 mag due to our 
Galaxy, obtained from the reddening map of Burstein & Heiles (1982). An additional 
LMC component of 0.11 mag was required, such that M v = —3.0 mag. In the absence of a 
far-UV extinction law, standard UV laws (Seaton 1979; Howarth 1983) are extrapolated for 
A < 1200 A with (variable) influence on the fit quality to FUSE data. 

Overall the observed spectrum of Sand 2 is very well reproduced by the model, with 
most Hell, Cm-lV and OlV-vi lines matched in strength and shape, except for Ovi 
AA3811-34 (model too weak) and Ovi AA1032-1038 (model too strong) as discussed above. 
In particular, the flat-topped nature of C ill A2297 is well matched. It is clear that although 
WO stars have little or no C ill A5696 (a classification diagnostic), other Cm lines are 
indeed present (Hillier 1989). Many spectral features in Sand 2 are due to blends because of 
the very broad spectral lines and the fact that recombination lines of O VI and C IV overlap 
with He II lines. For example, the spectral feature at AA4650-4686 has principal contributors 
Hell A4686, Cm AA4647-50, Civ A4658 and Civ A4685, while minor contributors include 
C IV A4646, A4689 and O vi A4678. 

Hell A5412/CIV A5471 provides an excellent diagnostic of C/He (e.g. Hillier & Miller 
1998) for spectroscopic studies of WC stars. However, the large line widths of WO stars, 
and the fact that Civ A5412 (14-8) contributes to Hell A5412 (7-4) hinders the use of 
these lines (see inset box in Fig. 2). Instead, we are able to derive C/He~ 0.7±0.2 by 
number from Hell (6-4) A6560/Civ A7700, although Civ (12-8) contributes to the former. 
Other Civ and Hell recombination lines show excellent agreement, except that Civ A1107 
is predicted to be stronger than FUSE observations reveal. 

In contrast to WC stars, numerous oxygen recombination lines are present in the 
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UV and optical spectra of WO stars (e.g. Ovi A5290, A2070). From their strength 
relative to He and C recombination lines, we estimate 0/He~ 0.151q;q5 by number. The 
weakness of lower ionization oxygen features, such as OlV A1400 and Ov A3150 also argue 
against higher O/He ratios, although Ov A5590 favours O/He~0.25. Poor fits to Ovi 
AA1032-38 and AA3811-34 are discussed above, while Ov AA2781-87 and Oiv AA3063-71 
are systematically too weak for all models. 

We present the predicted temperature structure of our Sand 2 model in Fig. 3. 
Although Sand 2 is an extremely hot star, its very high content of efficient C and O 
coolants (see Hillier 1989) directly results in a very cool (<10kK) outer wind (r/i?* > 100). 
Consequently, the ionization structure of metal species is predicted to be very stratified, 
such that 6+ is the dominant ionization stage for r/i?* <5, yet 3+ is dominant for 
r/R, >25. 

The high stellar temperature of Sand 2 implies a high bolometric correction (—5.7 
mag) and consequently hard ionizing spectrum, such that 50% of the emergent photons 
have energies greater than 13.6eV (912A Lyman edge), and 30% greater than 24.6eV (504A 
He I edge). The ionizing fluxes in the Hi, He I, and Hell continua are 10 49 1 s -1 , 10 48 8 s _1 
and 10 40 ' 3 s -1 , respectively. WO stars are known to produce strong nebular Hell A4686 
emission in associated Hn regions (e.g. Kingsburgh & Barlow 1995), although the predicted 
number of He II continuum ionizing photons in our Sand 2 model is much lower than those 
inferred from other WO stars. 



4. Discussion 



Abundances derived here qualitatively support the results from KBS, who used 
recombination line theory to derive C/He=0.5 and O/He=0.1 for Sand 2. WO stars are well 



- 9- 



suited to recombination studies for carbon, although oxygen is somewhat more problematic 
for recombination line studies, since the ionization structure is more complex (see Fig. 3), 
and few lines have available coefficients. 

Grafener, Hamann & Koesterke (1999) have carried out a detailed non-LTE 
spectroscopic analysis of Sand 2 (see Grafener et al. 1998). Overall, we confirm their M / yff 
determination, but derive a higher luminosity (by 0.2 dex) and temperature (they derived 
T*=101kK), attributable to the incorporation of line blanketing. More significantly, we 
obtain systematically lower metal abundances (they estimated C/He=1.3 and 0/He=1.2), 
such that the oxygen mass fraction for Sand 2 is only 16%, versus 50% according to 
Grafener et al. (1999). We attribute this major revision to improved spectroscopic and 
atomic datasets (Grafener et al. used simple C and O model atoms plus low S/N optical 
data). Our higher luminosity and clumpy wind conspire to revise the wind performance 
ratio, Mfoo/(L/c), from 56 to 12. 

Fig. 4 compares the luminosities and (C+0)/He ratios for Sand 2 with six LMC 
WC4 stars, updated from Dessart (1999). He improved upon similar work by Grafener 
et al. (1998) using line blanketed, clumped models, revealing a greater range of carbon 
abundances, 0.1<C/He<0.3 (due to improved spectroscopy), systematically higher 
luminosities (because of blanketing and improved reddenings), and lower mass- loss rates 
(due to clumping). The carbon enrichment of Sand 2 is substantially higher than the 
WC4 stars, although O/He does not differ so greatly from the WC4 sample, for which 
O/He<0.08. This supports the suggestion by Crowther (1999) that unusually high oxygen 
enrichment may not be a pre-requisite for a WO classification. 

We have superimposed (non-rotating) evolutionary tracks from Meynet et al. (1994) 
at 0AZ Q for M ini tiai=60, 85 and 12OM on Fig. 4. These evolutionary models predict 
C/0~2 when C/He~0.7, in conflict with our determination of C/0~4. Better agreement 
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is expected for evolutionary models in which rotation is accounted for, since these predict 
higher C/O ratios during the WC/WO phase (Maeder & Meynet 2000). From interior 
models, M initia i > 6OM for Sand 2, with a corresponding age of ~3-4.3 Myr, such that a 
supernova is expected within the next 0.1-5xl0 4 years. The stellar luminosity implies a 
current mass of 1OM (Schaerer & Maeder 1992), such that the mean post-main sequence 
mass-loss rate is 2.5-5xl0~ 5 M yr _1 . 

Quantitative analysis of other WO stars suffer from either (i) high interstellar reddening 
(WR102, WR142), (ii) complications because of binarity (Sand 1), or (iii) large distances 
(DR1). Nevertheless, we expect similar C and O enrichment to that derived here for Sand 2 
(KBS). 

This work is based on data obtained: (i) for the Guaranteed Time Team by the 
NASA-CNES-CSA FUSE mission, operated by the Johns Hopkins University, (ii) with 
the NASA/ESA Hubble Space Telescope obtained from the data archive at the STScI, 
which is operated by AURA under the NASA contract NAS5-26555; (iii) with the 2.3m 
Mount Stromlo and Siding Spring Observatory. Financial support is acknowledged from 
the Royal Society (PAC), NASA contract NAS5-32985 (U.S. FUSE participants), STScI 
grant GO-04550.01-92A, NASA grant NAG5-8211 (both JDH), PPARC (OD) and the UCL 
Perren Fund (LD). 
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Fig. 1.— Observed FUSE (CALFUSE 1.6.8) spectroscopy of Sand 2, rebinned to 15 kms" 1 , 
revealing the P Cygni VI resonance doublet at 1032-38A. The far-UV spectrum is affected 
by interstellar atomic and molecular hydrogen lines, plus airglow from [Oi], [Nl] and Ly/3. 

Fig. 2. — Comparison between far-UV (FUSE), UV (HST/FOS), optical and near-IR 
(MSSSO 2.3m) spectroscopy of Sand 2, shown as a solid line, together with our synthetic 
spectrum (dotted line), reddened by E(B - V)=0.08 (Gal) and 0.11 (LMC). Correction for 
interstellar Hi and H 2 in the far-UV has been applied following J.E. Herald, D.J. Hillier & 
R.E. Schulte-Ladbeck (in preparation), for which we adopt Tism=100K, t> tur b=10 kms -1 . We 
use log(n(Hl)/cm 2 )=20.6 (Gal) plus 21.3 (LMC) (Shull & van Steenberg 1985; Koornneef 
1982), plus log(n(H 2 )/cm 2 )=20. 

Fig. 3. — Theoretical temperature distribution for our Sand 2 model versus Rosseland optical 
depth, indicating selected stellar radii, i?*, plus the ionization balance of C, O and Fe. 

Fig. 4. — Comparison between the luminosity (in L Q ) and surface abundances ((C+0)/He, 
by number) of Sand 2 (triangle) with six LMC WC4 stars (circles, Dessart 1999), plus 
evolutionary predictions of Meynet et al. (1994) for M in i t i a i=60 (dotted), 85 (dot-dashed) 
and 120 M© (dashed), including locations of SN explosions (stars). 
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